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ABSTRACT 

We predict Lyman-a (Lya) luminosity functions (LFs) of Lya-selected galaxies (Lya 
emitters, or LAEs) at z = 3—6 using the phenomenological model of Dijkstra & Wyithe 
(2012). This model combines observed UV-LFs of Lyman-break galaxies (LBGs, or 
drop out galaxies), with constraints on their distribution of Lya line strengths as a 
function of UV-luminosity and redshift. Our analysis shows that while Lya LFs of 
LAEs are generally not Schechter functions, these provide a good description over the 
luminosity range of log 10 (L a /ergs -1 ) = 41 — 44. Motivated by this result, we predict 
Schechter function parameters at z = 3 — 6. Our analysis further shows that (i) the 
faint end slope of the Lya LF is steeper than that of the UV-LF of Lyman-break 
galaxies, (with a median otLya < —2.0 at z j>4), and [it) a turn-over in the Lya LF of 
LAEs at Lya luminosities 10 40 erg s -1 < L a <10 41 erg s -1 may signal a flattening of 
UV-LF of Lyman-break galaxies at —12 > Mjjy > —14. We discuss the implications 
of these results - which can be tested directly with upcoming surveys - for the Epoch 
of Reionization. 

Key words: galaxies: high-redshift - galaxies: luminosity function, mass function - 
cosmology: reionization ultraviolet: galaxies 


1 INTRODUCTION 

The luminosity function (LF) of galaxies provides one of the 
most basic statistical descriptions of a population of galax¬ 
ies. It describes the number density of galaxies in a given 
luminosity interval. Generally, the LF is well described by a 
|Schechter (1976) function 

L 

exp 


(p(L)dL = 




(i) 


with a normalization parameter </>*, an exponential cutoff 
at L > L*, and, a power law with faint-end-slope a for 
L L*. The parameters depend on wavelength considered, 
galaxy type (e.g., passive versus star forming), and cosmic 
time. 

At high redshift, galaxies are typically identified ei¬ 
ther through their broadband colors, for example using the 
drop-out or Lyman-break technique (Steidel et al. 1996), or 
through narrow-band searches aimed at detecting emission 


lines ( Partridge fe Peebles||196f Djorgovski et al.||1985 ) In 
particular, young star forming galaxies emit a significant 
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fraction of their radiation as Lyman-a (Lya) emission, and 
this method has been proved to be very efficient in finding 
samples out to z ~ 7 (e.g. 


Rhoads et al. 


2000 


Rhoads & 


Malhotra 2001; Ouclii et al. 2008 

Bond et al. 2009, 

2010 

Guaita et al. 2010; Kashikawa et a 

. 2011; Hibon et al. 

2012 

Ono et al. 2012 Ota & Iye|2012j R 

loads et al. 2012 ; Shibuya 

et al.|2012 Finkelstein et al.|2013 

Konno et al.||2014 

. 


Galaxies that have been selected (found) on the basis of 
their Lya lines are referred to as ‘Lya emitters’ (or LAEs). 
LAEs are useful because they are selected on having a strong 
Lya line flux irrespective of their associated UV-continuum 
emission. Therefore, LAEs can be fainter in the continuum 
compared to Lyman-break galaxies, and complement galaxy 
samples obtained via broadband searches which have been 
extensively carried out with the Hubble Space Telescope out 


to z ^ 

10 (e.g. Yan & Windhorst 2004 Beckwith et al. 

2006 

Bouwens^UilJ 2006 

Wilkins et al. 2010; Trenti et al. 

2011 

Grazian et al. 2012; 

finkelstein et al. 2012 

Bouwens et al. 

2014a 

|Finkeistein et al. 2014 Oesch et al. 

2014 Schmidt 

et al. 

2014). Moreover, the sensitivity of the observed Lya 


flux to intervening neutral hydrogen gas makes LAEs an ex¬ 
cellent probe of the Epoch of Rcionization (see e.g. [D ijkstra| 
2014 for a review). 
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Since the range of observed Lya luminosities at high -z 
typically extends only over ~ 1 — 1.5 orders of magnitude, 
the shape of the Lya LF is not strongly constrained and a 
fit with a Schechter function leads to significant degeneracy 
in the parameters. In particular the faint-end slope auyc, is 


essentially unconstrained: for example, Henry et al. (20121 
used a sample of six (three) LAEs to find a^ ya = — 1.70^'^ 
(QLyn = —1.45lo!?o) at 2 = 5.7. Other approaches include 
assuming a fixed value for OLya and resorting to the data 


to constrain the other parameters l 

van Breukelen et al. 

20051 Dawson et al. 2007, Ouchi et al 

2008 Hu et al. 2010 

Kashikawa et al. 2011, Ciardullo et al. 2012; Zheng et al. 

20131. In contrast, the ultraviolet (UV) LF of Lyman-break 


galaxies (LBGs) is much better constrained due to available 
data stretching over several orders of magnitude in lumi- 


nosity ( 

McLure et al. 2010 

Yan et al. 2011 

Bouwens et al. 

2011; Bradley et al. 2012 Oesch et al.|2012 

Yan et al. 2012; 

Lorenzoni et al. 

2013 

Schenker et al. 

20131. For the faint 


end slope, the most recent by Bouwens et al. (2014a| finds 
auv = —1.91 ±0.09 (auv = —1.64±0.04) at z ~ 6 (z ~ 4). 


There exists a clear opportunity to connect LAEs and 
LBGs via the Lya line emission properties of LBGs. |Shap-| 
ley et al. (|2003 1 provided a probability distribution func¬ 


tion (PDF) of the rest-frame equivalent-width (EW) of the 
Lya line in their sample of ~ 800 z ~ 3 LBGs. |Dijk-| 


stra & Wyithe (20121 showed that this observed PDF was 


well described by an exponential function, and that the 
characteristic scale-length of this function increased towards 
fainter UV-luminosities. While there do not exist equally 
well measured PDFs at higher redshifts and/or fainter UV- 
luminosities, recent studies have constrained both the red- 
shift and UV-luminosity dependence of the so-called ‘Lya 
fraction’, which quantifies the fraction of LBGs for which 
the Lya EW exceeds a certain value. The Lya fractions - 
which represent integrated versions of the full EW PDF - 


2010 Cassata et al. 2015) and from UV-bright to UV-faint 

galaxies (Stark et al.|2010 12011 

Pentericci et al. 2011 Ono 

|et al. |2012, Schenker et al. 2012 

)• 


There have been several attempts to link the redshift 
evolution of LBGs and their Lya fractions to LAE luminos¬ 


ity functions (Dijkstra & Wyithe 2012 Faisst et al. 2014 


Schenker et al. 20141. In this paper, we follow the work 


of Dijkstra & Wyithe (20121 and combine the most recent 


constraints on UV-LFs & Lya fractions to make predic¬ 


tions for Lya LFs. Dijkstra & Wyithe (20121 showed that 


this phenomenological model reproduces observed Lya LFs 
and their redshift evolution remarkably well. Here, we focus 
specifically on the faint-end slope of the Lya LF of LAEs, 
because (i) we can make robust predictions for this faint end 
slope, (it) as we will show later, this faint end slope can be 
highly relevant for understanding the Epoch of Reionization. 


This paper is structured as follows. In Sec. [2] we lay out 
our method. We present our results in Sec. IsTand discuss 
them in Sec. [4] Finally, we conclude in Sec. [5j The cosmo¬ 
logical parameters we adopt are D m = 0.3, Da = 0.7, h = 
0.7, a 8 = 0.9. 


2 METHOD 

The number density of LAEs with luminosities in the inter¬ 
val [L a ± dL Q /2] is given by 

M V V, max 

4>hAE(L a )dL a = dL a F j dMuv 4 > (Muv)P(La\Muv) 

MuV, min 

( 2 ) 

Here, <j>(Muv)dMuv denotes the number density of LBGs as 
a function of in the range Muv±dMuv/2. This function can 
be represented by the Schechter function with parameters 
(auv , M^jy, (j>uv)- 

The term P(L a \Muv)dL a is the conditional probabil¬ 
ity that a galaxy has a Lya luminosity L a given an abso¬ 
lute UV magnitude Muv- This conditional probability can 
be recast in terms of the equivalent width (EW) prob¬ 
ability density function P(EW\Muv) as P(L a \Muy) = 
P(EW\M UV )^^ if EW > EWuae, where L a and EW are 
related as L a = EWL\ = EW[vL v / A]. Here, the luminos¬ 
ity/flux densities, frequency and wavelength are evaluated 
just longward of the Lya resonance at A = (1216 ± e)A. We 
can extrapolate these flux/luminosities to their values where 
the UV-continuum measurements are usually made (see e.g. 
Dijkstra & Westra 2010)Q Furthermore, EWlae denotes 
the equivalent width threshold that determines whether a 
galaxy would make it into an LAE sample. We adopt that 
EWlae = 0A, but note that some surveys adopt colour 
criteria for selecting LAEs as large as EITlae = 64A 
(see Dijkstra & Wyithe 2012). If EW ^ -EWlae, then 


P(L a \Muv) = 0 since in this case the galaxy does not qual¬ 
ify as an LAE. This threshold more closely represents de¬ 
tection threshold for Lya emi tting galaxies in sp ectroscopic 
survey^]- e.g., with MUSE ( Bacon et al ||2010|), HE TDEX 
( |Hill et al.|2008 1 and/or VIMOS ( Cassata et al.|2011 20151. 
We have verified that our main results do not depend on this 
choic<j3 

The preceding factor F in Eq. is merely a normaliza¬ 
tion constant to fit the data and, hence, can be thought of as 
the ratio of predicted versus the total number of LAEs. This 
factor should ideally be F = 1. However, Dijkstra <V W yithe| 
(20121) required that F ~ 0.5. The origin of this number is 


not known (see Dijkstra & Wyithe 2012 for an extensive 
discussion J] but we stress it only affects the predicted nor¬ 
malization linearly and not the predicted faint-end slopes. 
Hence, the key function in our analysis is P(EW\Muv)- 


1 We use the relation La. = C\EW Ljjy u where Ljjy u oc 
v~P~ 2 is the UV luminosity density Ljjy u and C\ = 
Va/\a(Xuv/X a )~^~ 2 converts the flux density at A = (1216 + 
e)A to that at Xjjv = 1600A, which is the wavelength where 


Ljjv,v was measured (Dijkstra & Westra|2010 I. 

2 Note that in practise an EW cut is likely still needed to distin¬ 
guish between LAEs and lower- 2 : interlopers, such as [Oil] emit¬ 
ters. This EW cut can nevertheless be lower than EWlae ~ 20A 
(Leung et al.|2015) 

^ We have verified that varying EWlae in the range [0, 50] A 
changes QfLyo: by ~ 0.02. 

4 The value of F depends weakly on the adopted UV Schechter 
function parameters. For example, |Bowler et al.| ( |2014[ ) reported 
slightly different best-fit values, which drive F up to F ~ 0.7—0.8. 
The |Finkelstein et al.| ( |2014| ) parameters, on the other hand, also 
suggest F ~ 0.5. 
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Several functional forms have been explored in the literature. 


Schenker et al. (2014) compared the maximum likelihood 


values for several EW distributions to their Keck MOSFIRE 
( McLean et al.|2012 1 data, and concluded that the exponen¬ 
tial distribution introduced by Dijkstra & Wyithe (20121 


provides an adequate fit. This functional form is 


P{EW\M uv ,z) = A/" exp 


EW 


EW C {M UV , z) 


(3) 


with EW C — E]V C , o d - yM uv (Muv T Muv, o) T Pz(,z + zo ) 
where Pm uv , Pz, Muv, o, 2o, and EW C , o are model param¬ 
eters. These parameters were chosen to match the observa¬ 
tions of ( Shapley et al.|2003 1 and Stark et al. (2010 2011 \ as 
closely as possible. Furthermore, AT is a normalization con¬ 
stant which is forced to be zero outside of [EW m i n , EWmax]- 
Our choice of values for the model parameters is described 
in Sec. 3.3| where we present the numerical results. In Ap¬ 
pendix Slwe show explicitly that the main results in this 
paper are insensitive to both the functional form of P(EW) 
and the parameterization of EW C - 


3 RESULTS 


We first present results in which EW C =constant (in § 3.1l. 


This allows us to demonstrate that for models in which the 
Lya fraction does not evolve with Muv, the faint end slope 
of the LF of LAEs approaches that of LBGs. We then present 
a simplified model in § |3.2| in which the mean Lya EW- 
PDF increases towards fainter UV-luminosity function. This 
model demonstrates quantitatively that the faint end slope 
of the LF of LAEs is steeper than that of LBGs if the Lya 
fraction increases towards fainter UV-luminosities. In § |3.3| 
we present the results that we obtained from the EW-PDF 
given in Eq. 


3.1 Exemplary case with EW C = const. 

We consider the case P(EW\Muv ) = P{EW), i.e. EW C = 
const. = X/Ci. Furthermore, we set Af = 0 for EW < 
EWlae- Under these assumptions we find 


OO 



where K n {x) is the modified Bessel function of the second 
kind and L^ v is the luminosity corresponding to M^ v . 

Eq. § shows that the Lya LF generally does not take- 
on a Schechter form. The slope of the LF is given by 

_ d log 4>{L a ) _ FyK a UV— 1 Qy/V) (a 

“ Ly “- d log L a K auv (2^/y) U 

with y = L a /(Lu V A). For L a Lu V \ we have y « 1, and 
we obtain to leading order aLya « — r(l— auv)/F(— oijjv) = 
auv- Thus, having a constant EW C corresponds to an un¬ 
changed faint end slope, auya = auv- 



Figure 1. Upper panel: The predicted number of LAEs in the 
range log 10 (L a /er gs" 1 ) ± dlog 10 (Lq/erg s -1 ) using the UV 
LF evolution from |Bouwens et al.] ( |2014a| ) taken at z = 5.7 
(black solid line). The grey dashed line marks the faint end slope 
(aLya — —1.90) and the dashed dotted lines show Schechter 
fits to our numerical findings. Once the fit was carried out 
over the whole shown luminosity range (blue) and once only 
in log 10 (L/ergs -x ) = [42,43.5] (green). The red discs are the 
z = 5.7 observations by |Ouchi et al. C2 008[ > and the black arrows 
denote the MUSE deep field and medium deep field as well as the 
JWST limits at that redshift (see text for details). Lower panel: 
Relative deviation of the fits to the numerical results. 


3.2 Exemplary case where P(EW\Muv) evolves 
with Muv 

If EW C depends on Muv a general analytic solution for 
<f>(L a ) does not exist. For illustration purposes we first con¬ 
sider a case in which we replace Eq. (J3]) with a Dirac-5 dis¬ 
tribution, 

p(EW\Luv) oc S(EW - EWd) , (7) 

where EWd(Luv) = EWd,o {Luv / LuvE■ The parameter 
EWd can be interpreted as the mean of the full PDF. This 
5-function PDF lead0to 


4>{L a )AL a oc AL a x L“ j/v/(7+1) 

_ ( L a 


x exp 


KCxEWd'OLt 


1 /( 7 +!)' 


• ( 8 ) 


Here, the faint end slope is aL yQ = ctuv /{7 + 1)- he Lya 
LF thus has a steeper faint-end slope than the LBG LF, 
if 7 < 0 (i.e. if EWd decreases towards fainter Luv, as has 
been observed). Also note that we again obtain auya = oiuv 
if EWd does not evolve with Muv- 


3.3 Realistic case with P(EW\Muv) inferred from 
observations 


For the model parameters of P(E!F|Muy) in Eq. we 
adopt the values from Dijkstra & Wyithe ( 2012[ p| Example 


5 For simplicity, we set the minimum and maximum UV lumi¬ 
nosity to zero and infinity, respectively. 

6 Specifically the model parameters related to EW C 

are given by (EW C)0 , ii Muv , Mz, M UV)0 , z 0 , F) = 
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Figure 2. Evolution of the Schechter parameters of the UV LF from |Bouwens et al.| ( [2014a| ) (as blue dashed lines) and the parameters of 
the computed Lyct LF ( black lines). In particular, the left panel shows the characteristic luminosity L* (note the different normalization 
constants), the central pane 1 the faint end slope a, and the right panel the overall normalization (j>* . Predictions at z > 6 (within the shaded 
grey area) do not account for reionization (see § |4.3| ) ■ In this region, the black solid lines correspond to models with an uninterrupted 
EW evolution, whereas the dashed-dotted lines represent a model in which we freeze the EW evolution, i.e., EW c (z > 6) = KW r {z = 6) 
(this assumption has been adopted in previous works). 


EW-PDFs are shown in Appendix |A1| For a more detailed 


motivation of this P(E\Y) we refer the reader to Dijkstra 


& Wyithe (2012). We integrate the UV-LF over the range 


Muv,( min,max) = (—30, —12) when predicting Lya luminos¬ 
ity functions, and discuss the impact of varying M max in 
Sec. |4] 

The redshift evolution of the best fit Schechter param¬ 
eters of the UV LF is taken from Bouwens et al. (2014a) 
and given as Mf, v = —20.89 + 0.12 z, <j>uv = 0-48 x 
10 -O.19(*-6) 10 —3 cMpc — 3 J andj auv = _ L85 _ 0.09(z - 6). 

Following these analyses, we use Xuv = 1600 A as rest frame 
wavelength in which the UV continuum was measured and 
assume a UV spectral slope /3 = —1.7. This choice for /3 
does not affect our results (see Appendix A4 for detailed 
discussion). 

The upper panel of Fig. [T] shows the resulting number 
density of LAEs at z = 5.7 in the luminosity range log 10 L a ± 
dlog 10 L a /2, i.e., ip(L a )d log 10 L a , as a function of L a . This 
quantity is related to 4>(L a ) as ip(L a ) = cf){L a )L a log 10 
(‘log’ denotes the natural logarithm). We compare these pre¬ 
diction to the data from Ouchi et ah (2008). In addition, we 
show the MUSE detection limit^for its medium deep field 
(MDF, limiting flux F > 1.1 x 10 -18 ergs -1 cm -2 , integra¬ 
tion time Tint. = 10h), and, deep field (DF, F > 3.9 x 
10 -19 ergs _1 cm -2 , Tint. = 80h) surveys as well as an exem¬ 
plary JWS"#]limit (F > 10 -18 erg s -1 cm -2 , T int . = 10 4 s). 

Figure. [Halso shows two Schechter function approxima¬ 
tions to our numerical findings fitted over the full luminosity- 
range shown (in blue) and over log 10 (L a /erg s" 1 ) = 

[40.5,42] (in green). Although we do not expect the resulting 
Lya LF to be a Schechter function (as shown in Sec. |3.2[ ), it 
provides a reasonable fit over the displayed luminosity range. 


(23A,7A,6A, 21.9, -4.0,0.53). The EW-PDF covers the 
range [E W m i n , Elfmax] ■ Here, the lower limit EW m j n = —ai, 
where ai(Muy) follows the form a \ = 20 A for Mjjy < —21.5, 
ai = (20 - 6 (M uv + 21.5) 2 ) A for -21.5 sf M uv sf -19.0 and 
ai = —17.5 A, otherwise (see Dijkstra & Wyithe 2012). We used 
BWmax = 1000 A but we verified that this choice does not affect 
our results quantitatively. 

7 MUSE survey limits taken from http://muse.unlv-lyonl.fr/ 
IMG/pdf/science_case_gal_formation.pdf 

° JWST survey limits obtained from http://www.stsci.edu/ 
jwst/science/sensitivity/ 



- 39 

— 40 

— 41 

— 42 


Figure 3. p(Muv\L a ) oc p(L a \Muv)<t>{Muv) versus Muv for 
some exemplary values of L a . The black arrow shows the JWST 
photometric limit quoted by |Windhorst et al. | |2006| for z ~ 6 
which corresponds to T; nt . ~ 10° s integration time. 


This can also be seen in the lower panel of Fig. [l] where we 
display the relative deviation of the fits to the LF. 

Fig.[2]shows the redshift evolution of the Schechter best 
fit parameters (as black lines). Predictions for z > 6 do not 
account for reionization effects and are calculated with an 
unaltered EW evolution (solid line) as well as an EW-PDF 
which does not evolve after 2 = 6 (dash-dotted line). We 
discuss this result separately in § |4.3| For comparison, we 
plot the corresponding redshift parameterization of the UV 
LF by Bouwens et al. (2014a I (as blue dashed lines). The 
left panel shows that L* increases by a factor ~ 2 over the 
redshift range 2 = 3 — 6, which differs from the redshift 
evolution in the characteristic UV-luminosity which drops 
by 20%. This difference is driven by the redshift evolu¬ 
tion in the Lya EW-PDF, which in turn was inferred from 
the observed redshift-evolution of Lya ‘fractions’ over this 
redshift range. The central panel shows that aLy a < auv. 
This is again a consequence of inferred redshift evolution 
of the Lya-EW PDF (see § 3.21. This figure also illustrates 
the close-to-linear ay ya -z relation. This evolution is mostly 
driven by the redshift evolution of auv- Finally, the right 
panel shows the predicted redshift evolution in <f>*. 
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4 DISCUSSION 
4.1 Low -L turnover 

The integral over (j>{L a )AL a diverges for a < —2. We there¬ 
fore expect that the luminosity function flattens or turns- 
over below some luminosity. The minimum luminosity that 
we can account for in our models is 


= £W m in(MuV ,max ) CiLuv ,n 


(9) 


where -EWmm = — ai = 17.5 A (see footnote [ii] in § 3.31 


denotes the minimum equivalent width in our EW-PDF at 
the maximum absolute UV-magnitude (i.e. the lowest UV- 
luminosity). For example, we obtain L ajm i n ~ 10 39 ergs _1 
for Muv.max = —12. At this luminosity we expect the pre¬ 
dicted Lya luminosity to go to zero, as is shown in Figure]!] 

An estimate for where we may start to see departures 
from a power-law slope can be obtained by considering the 
conditional probability p(Muv\L a ). Bayes’ theorem states 
that p(Muv\L a ) oc <f>(L a \Muv)(f>(Muv ), of which we show 
examples in Fig. [3] for four different values of L a . This Fig¬ 
ure illustrates for example that Lycr observations that probe 
a flux corresponding to L a = 10 4 ° erg s _1 - a level that can 
be reached in MUSE ultra deep fields - effectively probe 
galaxies with —14 < Muv < —11, which are fainter than 
can be probed directly even with the JWST. The JWST 
detection limit shown in Figure [3] is taken from |Windhorst| 
|et al.| ( |2006| ). Figure [3] further shows that if the UV-LF flat¬ 
tens off at - say - Muv )> — 12 that then the effects should 
become noticeable in the predicted Lya luminosity function 
around L a = 10 4t) erg s _1 , as here galaxies with Muv ~ —12 
dominate the contribution to the Lya LF. 

In Figure [4] we make these points more explicit, and 
show the predicted faint end of the LAE LF for four val¬ 
ues of Muv,max (calculated with the UV LF parameters at 
z = 3.1). For each curve we marked L a<m in with dotted lines. 
For example, a potential UV turnover at Muv ~ —12 leads 
to deviations]^] of the Lya LF at L a ~ 10 4 ° ergs -1 and a cut- 

also contains data points 


off at L a 


10 


l39 erg s 1 . Figure |4 


taken from [Rauch et ah (2008). 


Rauch et al. (2008) per¬ 


formed an ultra deep (92-hr) exposure with VLT s FORS2 
low resolution spectrograph. The goal of these observations 
was to detect fluorescent Lya emission from optically thick 
clouds powered by the ionizing background. While their sen¬ 
sitivity turned out not to be good enough to detect this 
fluorescent emission (revised estimates of the ionizing back¬ 
ground and the conversion efficiency into Lya), they de¬ 
tected numerous ultra faint Lya emitting sources charac¬ 
terizing their LF down to L a ~ 6 x 10 4O ergs -1 . We com¬ 
puted the uncertainties with the cosmic variance calcula¬ 


tor of Trenti & Stiavelli (20081. These data-points fall on 


the predicted LF for Muv.max = —16. However we caution 
that the turn-over occurs at the lowest luminosity data-point 
only, which might suffer from incompleteness (although it 
lies above the detection threshold). In the same figure we 
provide the estimated MUSE limits for the deep field (DF) 
and the gravitationally lensed ultra deep field (UDF) sur¬ 
veys. 



Figure 4. ip(L a )dL a versus L a at z = 3.1 for different values 
of Muv max illustrating the cutoff at low Lya luminosities. The 
dashed and dotted lines show the “cutoff” and “deviation” points 
for each Muv ,max discussed in Sec.[4| Data points are from|Rauch| 
|et al.| (j2008| taken at that redshift(see the text for a discussion 
on the error bars). And the black (grey) arrow denotes planned 
future MUSE (ultra) deep field limits. 


4.2 Implications for the Epoch of Reionization 


Low luminosity galaxies are expected to play a major role in 
driving the reionization of the Universe (e.g. Robertson et al. 
2010 ; |Trenti et al.|[2 01 0] |Kuhlen fe Faucher-Giguere||2012| 


Boylan-Kolchin et al. 2014). Determining the faint end of the 
LF such as its slope and a turnover luminosity is essential 
for constraining the volume emissivity of ionizing photons. 
However, even future experiments will have difficulties de¬ 
tecting these galaxies directly via their UV continuum flux. 
Current constraints rely, therefore, on extrapolation of local 
properties to higher redshifts (Weisz et al. 2 014]) , (relatively 
few) gravitationally lensed objects ( Alavi et al.||2014 |Atek 
et al.|20l4 I or inferences from gamma-ray burst observations 
( |Trenti et al. 20121. In this work, we have shown that the 
Lya LF can provide an independent probe of the faint end 
of the UV LF, and that for example the MUSE DF survey 
could already detect (or rule out) a turnover at Muv —15. 

Recent studies have shown that Lya escape may be cor¬ 


related with the escape of ionizing photons (Behrens et al. 
2014 Verhamme et al. 2014), as the escape of ionizing pho¬ 
tons requires low Hl-column density (Nm < 10 17 cm -2 ) 


channels, which can also provide escape routes for Lya pho¬ 
tons. The fact that Lya LFs are likely steeper than the UV 
LFs implies that the Lya volume emissivity - and therefore 
possibly the ionizing emissivity - are weighted more strongly 
towards low luminosity galaxies. This is consistent with the 
expectation that ionizing photons escape more easily from 
lower mass - and hence lower luminosity - galaxies. A steep 
faint-end slope of the Lya LF may therefore provide obser¬ 
vational support for this scenario. 


4.3 Predictions for redshifts 2 = 6 — 8 

We extrapolated our predictions for the best-fit Schechter 
parameters of the LAE LF to 2 > 6 in two ways (shown in 
Fig# (%) in the first, we assume that the EW-PDF contin¬ 
ues to evolve as inferred from the observations at z = 3 — 6. 


9 The first deviations in the Lyo LF can be found at L a dev. = 
CiLuy(x)EW c (x) with x = Mjjy , m ax — A Muy. Here, A M\jy 


describes the half width of p(Mjjy\L a ) at a chosen probability 
threshold. 
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This model is represented by the solid lines, and, (ii) in the 
second, we ‘freeze’ the EW distribution for z > 6 at the 
value it had at a = 6 (dashed lines). This latter assumption 


has been common in previous works (see e.g. 

Dijkstra et al. 

201 lj Bolton & Haehnelt 2013; Jensen et al. 

2013~[Choud- 

hury et al.j!2014 

Mesinger et al. 12015 

). We show results for 


these two models to get a sense for the uncertainties on our 
predictions. We stress that we have purposefully not mod¬ 
elled the impact of reionization on the EW-PDF. Reion¬ 
ization is likely responsible for the observed ‘drop’ in the 


observed Lya fractions at z > 6 (e.g. Pentericci et al. 

|2011 

Schenker et al.|2012; Ono et al. 2012||Treu et al. 2013 

Caru- 

ana et al.|2014| |Tilvi et al. 2014). Understanding this drop 


has been the main focus of previous works, and is outside 
the scope of this paper. Our predictions for z = 6 — 8 are 
useful in a different way, as they provide predictions for the 
Lya LFs of LAEs in the absence of reionization. Compari¬ 
son to observed LFs at these redshifts highlight the impact 
of reionization. 


5 CONCLUSIONS 

We predicted Lya luminosity functions (LFs) of Lya- 
selected galaxies (Lya emitters, or LAEs) at z = 3 — 6 
using the phenomenological model of |Dijkstra fe Wyithe[ 
([2012). This model combines observed UV-LFs of Lyman- 
break galaxies (LBGs), with observational constraints on the 
Lya EW PDF of these LBGs, as a function of Muv and red- 
shift. The results from our analysis can be summarized as 
follows: 


• While Lya luminosity functions of LAEs are generally 
not Schechter functions, these provide a good description 
over the luminosity range of log 10 (L a /ergs -1 ) = 41 — 44 
(see Fig. |T|) . 

• We predict Schechter function parameters at a = 3 — 6 
(shown in Fig. [2|. The faint end slope of the Lya LF is 
steeper than that of the UV-LF of LBGs, with a median 
ctLyo, < —2.0 at z )>4 (see the central panel in Fig. 0. While 
the current work was in the advanced stage of completion, 


Dressier et al. (2014) posted a preprint in which they ob- 
servationally infer a very steep faint end slope at z ~ 5.7 
(—2.35 < a < —1.95, also see Dressier et al. (2011)). The 
central value a = —2.15 is in excellent agreement with the 
value a ~ —2.1 predicted in our framework. 

• The faint end of the LAE LF provides independent con¬ 
straints on the very faint end of the UV-LF of LBGs. For 
example, the predicted LAE LF at Lya luminosities 10 40 erg 
s -1 < L a <10 41 erg s -1 is sensitive to the UV-LF of LBGs 
in the range —11 > Muv > —15 (see Fig. [3] and Fig. [4j). 
These LBGs are too faint to be detected directly (even with 
JWST). A turn-over in the Lya LF of LAEs may signal a 
flattening of UV-LF of LBGs. We discuss implications of 
these results for the Epoch of Reionization in § |4.2| 


We have verified that these results are insensitive to our 
assumed functional form of P(EW) and how we parameter¬ 
ized its dependence on z and Muv- Our predictions can be 
tested directly with various upcoming surveys. 
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APPENDIX A: VARYING THE EW 
DISTRIBUTION 

In this appendix, we demonstrated that our main results 
and conclusions do not depend on our assumed EW-PDF. 


Al Fiducial P(EW) 

Our default EW distribution is given by Eq. @. We plot the 
PDF for three UV magnitudes and two redshifts in Fig. |A1| 
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A2 Schenker et alT] ( |2014 ) parameterization 

As mentioned in Sec. [2| [Schen ker et al. (20141 suggest an 
alternative parameterization of the EW-PDF, namely 


P(EW\/3) = 




: exp 


(log Ew-^m 2 


2 a 2 


(Al) 

a and 


V2naEW 

This log-normal PDF possesses the parameters A em , 

/r. The latter is given by 

/z(/3) =/xa + - 2.0), (A2) 


where j3 is the UV continuum slope. Schenker et al. (2014) 


found their EW distribution to depend more strongly on /3 
than on and therefore constrained P(EW\/3). We 

can include this parameterization into our formalism if we 
map P(EW\/3) onto P(EW\M vv ,z). 

This mapping is based on three results from Bouwens| 


et al. (2014b): 


(i) We use their empirical linear correlation between /3 — 
Muv at 2 ~ 7. This relation constrains /3(Muv = —19.5) = 
-2.05 ±0.09 ±0.13. 

(ii) Furthermore, we apply their measured change per 
unit redshift A/3 /Az = —0.1 ± 0.05. 

(iii) Finally, we use their measurement that A/3/AMuv = 
-0.2 (-0.08) for Muv -19 (> -19). 

Accordingly, our mapping can be written as 

0 = Po ± (M uv ± 19) + (A^ {z — 7), (A3) 

where Am vv = A/3/AMuv and / A^ = A/3/Az. 

In Eq. we used the best parameters by 

Schenker et ah] (20141, i.e., (A cm , a, fj, a , fj, s ) = 
(1.0, 1.3,2.875, —1.125). In addition, since A cm is de¬ 
generate with F, we set F — 1. The orange dashed line 
in Fig. |A2| shows the resulting Lya LF at 2 = 5.7. The 
agreement in the faint-end between the two procedures is 
remarkable. For greater luminosities, however, the (Schenker] 
et al. (2014) parameterization leads to a (much) higher 
number density of LAEs. While there are significant 
uncertainties in the above procedure, the agreement we 
get at the faint end slope is especially encouraging. Future 
surveys can be extremely useful in further connecting the 
LAE and LBG populations by constraining the bright end 
of the LAE luminosity function. 


A3 ‘Freezing’ the EW PDF for faint galaxies 

Since the evolution in the EW PDF for fainter sources in¬ 
volves a (modest) extrapolation of observationally inferred 
P(EW), we have also tested and alternative PDF where 
we ‘froze’ the evolution at Muv = —19. That is, we also 
conservatively assume that the EW-PDF stops evolving at 
Muv > —19 (even though observations hint that this is not 
the case, see fig. 13 of |Stark et al.| ( [201 0|). Fig. |A2| shows 
the resulting Lya? LF (purple line). It is clear that our re¬ 
sults are only affected slightly, i.e., the faint-end-slope «Ly Q 
is reduced by ~ 0.05. We also tested this over a variety of 
redshifts. 



Figure A2. The Lya LF at 2 = 5.7 for different EW distribu¬ 
tions. The solid black line and Ouchi et al.| |2008j data points 
are the same as shown in Fig. [T] As com parison we show the Lya 
LF computed using the |Schenker et al. | (|2014) parameterization of 
the EW PDF once with the full /3 —Muv relation as given in § |A2| 
(blue dashed line) and once with a constant Am vv = ~ over 
all Muv (green dashed-dotted line). Also shown are our results 
when freezing the EW PDF at Muv = —19 (purple) and when 
using the /3 — Muv relation instead of a constant 6 (orange). 


troduces some additional dispersion in the predicted Lya 
flux at a fixed Muv- However, varying /3 within [—2.0, —1.5] 
changes the Lya flux only by 1 — (Auv/ALya ) 1 ' 5-2 0 ~ 13%. 
This dispersion is smaller than that introduced by the EW- 
PDF. If we replace the constant /3 with the empirical fit 
described in § A2 then our predicted Lya LF (represented 
by the orange line in Fig. A2 ) is barely any different from 
our fiducial model that used /3 = —1.7 (represented by the 
black solid line). 


A4 Non-constant UV spectral slope 

The spectral slope /3 is not a constant, but depends on 
UV magnitude and redshift (as discussed above). This in- 
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